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Nomenclature 

A Area (m 2 ) 

C p Specific Heat Capacity, for Air (1005 J/Kg.K) 

K Thermal Conductivity (W/m.K) 

m Mass Flow Rate (Kg/sec) 

P Pressure (Pa) 

Pr Prandtl No., for Air (0.7) 

q Heat Added (J/Kg) 

R Gas Constant, for Air (287. 1 J/Kg.K) 

T Temperature (K°) 

V Velocity (m/sec) 

p Density (Kg/m 3 ) 

[i Dynamic Viscosity (Kg/m. sec) 
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Abstract 

A new method has been proposed for conversion the wind energy into a 
kinetic energy using the hybrid effect of nozzle and solar energy. As the air 
flows through the horizontal nozzle it accelerates because of the narrowing 
constriction. Heat energy acquired from the wall is converted into the kinetic 
energy of flow. Critical dimensions are calculated for the convergent nozzle 
that made of steel. This study focuses up on the benefits of using solar 
nozzle on the wind energy. The study is an attempt to raise the local wind 
velocity (3 m/sec) to a high velocity gives good energy allowed the wind 
plant to generate power, on other words in order to increase the efficiency of 
the plant. From the obtained results, it is observed that the velocity of wind 
increased by the increment of heat gain and decrement of the area. The 
velocity value in the case of heat added is reached to (50 m/sec), while in the 
case of no heat transferred is about (48 m/sec). Calculation indicates that 
maximum heat gained could give (57 KW) output power. 



J ^Imtj jliium AjS^ya. <iUa ^Lj^ll <SUa Jjj^J ^3 ^ij^ 3 jj;^-" Cllaall I iA ^ 

Afr^uill <Jx2 ^Lj^l ejLjj ^ -C5^ J^-"-^ ^ y-^ 3 J* ls^-"^ ( -5 -'-"-n''' j» J^ 3 ^- 

<c JjSj Uiu ^ jiJI ^oc JU. ^ (48 m/sec) *l £ jj=* ^ uj^ (m/sec 
Aiuiaj 3JU ^ .(jJa^ll J^i\ ^j^^ill ^ jJa j^JI JU. ^ (50 m/sec) £JJ=* 

.(57 KW) U jl^Lo ajS 4iLL ^ Jj.^ti ijUV^ti 'uUc 
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1. Introduction 

The improvement of any power plant needs to know what the components 
that depending on it. In the case of wind plant, it is necessary to know the 
behavior of the wind. There are several factors affect on the wind state such 
as velocity, temperature and flow direction. Power of the wind can be 
extracted by allowing it to blow past moving wings that exert torque or 
rotor. The wind profiles show that the wind energy conversion device with a 
certain amount of confidence to extrapolate better wind speed taken at a 
height more than (10 m) [1]. 

Solar power is the technology of obtaining usable energy from sun light. The 
advantages of used renewable energies (solar and wind) are; no pollution, no 
fuel needed and low costs [2]. It is clear that use of the solar energy at 
latitudes with higher levels of radiation will improve the performance of 
wind system and make it possible to cover loads with smaller and less 
expensive systems. Attempts to parameterize the solar radiation received at 
the ground are made for the global radiation (direct plus scattered, and for 
tilted planes also radiation reflected onto the surface), rather than separately 
for normal incidence and scattered radiation [3]. 

2. Background 

Padki and Sherif 1999 [4] suggested a small model solar chimney of 
radically different geometry. In this case the chimney is conical with the 
turbine placed at the narrowing of the cone. The paper is mainly theoretical 
but they argue that there is no upper bound on the efficiency. The efficiency 
is rise up to 20% by the heat gained. Unfortunately the practical results 
reported are much less impressive. 

An international study of power system impacts of wind power [5] has been 
formed in 2006 under the IE A implementing agreement on wind energy; the 
task will analyze case studies from different power systems. The results are 
not easy to compare, for example the incremental regulation due to wind was 
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found to be 36 MW in a location in USA while in another location in 
Netherlands it is about 6000 MW. 

Williams 2007 [6] proposes theoretically a prototype convergent nozzle 
made of glass and arranged vertically. A solar absorber abundantly 
perforated and several layers thick or of metallic honeycomb structure is 
arranged in the lower levels but above the base of the nozzle. The 
convergent solar nozzle has base diameter of (10 m), throat diameter of (1 
m), and height of (10 m). Calculation indicates that the maximum output 
power in the summer of UK is about (59 KW) and the temperature 
difference through the nozzle is (1.2 °C). 

Sakonidou and et al. (2008) [7] determine the tilt that maximizes natural air 
flow inside a solar chimney. The model starts by calculating the solar 
irradiation components absorbed by the solar chimney. The calculations 
have been obtained by a simulation program solves the relevant conservation 
of mass, momentum and energy equations (CFD) for solar chimney using 
finite difference methods. The model predicts the temperature and velocity 
of the air inside the chimney as well as the temperatures of the glazing and 
the black painted absorber. The experimental chimney duct has the shape of 
a narrow parallelepiped with dimensions: 1 m height, 0.74 m width and 0. 1 1 
m gap. Black painted aluminum sheet is used for the construction of walls of 
the chimney. The walls have absorptance of (0.95). Comparisons of the 
experimental model predictions with CFD calculations delineate that there is 
a good agreement between them at different tilt positions. 

The objective of the present work is to present a theoretical analysis for the 
using of solar nozzle. The solution is merged between analysis of wind 
energy and solar energy to increase the wind energy to a rate allowed the 
wind plant to generate power, also the efficiency of the plant can be 
increased by using a convergent nozzle arranged horizontally and exposed to 
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solar radiation. A wind turbine placed in the throat of the nozzle converts 
flow kinetic energy into electricity, see figure (1). 



3. Solar Radiation Modeling 

The fact that the amount of solar radiation received at a given location at the 
Earth surface varies with time due the Earth's rotation (Diurnal Cycle) and 
also depending on the latitude [3]. The average rate of the solar radiation 
incidence on a horizontal plate in Iraq for July is about (450 W/m 2 ) [1]. 

Suppose that the nozzle we had is the solar collector, so the input radiation 
power is [8]: 

q R =I N ocx (1) 

I N =I s Cos9 (2) 

Cos 9 = Cos 4» Cos 8 + Sin § Sin 8 (3) 

Where: 

ax = Product of absorptive and transitivity (0.95 [7]) 

Is = Solar Radiation (450 W/m 2 ) [1] 
For Baghdad (<)> = 33°) at noon time in July (8 = 23.4° ) [3] 
The incident angle is nearly vertical so (cos 6 = 0.97). 



q 




Figure (1) Supposed Wind Plant 
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There is heat loss across the wall denoted as [8]: 

q L =U(T b -T a ) —(4) 

Where T a = Air temperature (35 °C) 

T b = Average bulk temperature inside the nozzle (about 36 °C) 

Principally, the overall heat transfer coefficient (U) can be calculated from 
[9]: 

U = T ^ --(5) 

Where, i:- refers to the internal part of the nozzle 
o:- refers to the external part of the nozzle 
The convection heat transfer coefficient (h) can be determined from [9] : 

h = ^ --(6) 
d 

Nusselt No. and Reynolds No. can be determined from [9]: 
Nu = 0.023 Re 08 Pr 03 (7) 

Re=P^ --(8) 

Consequently, the useful heat gained to the flow will be: 

q = a n f r (q R -qj ( 9 ) 

The heat removed factor (F R ) has value of (0.8) [8]. 

A N is the surface area of the nozzle which represents a semi-cone shape so it 
is calculated from: 

A N =rc^±^L ---(10) 

Where, L is the nozzle length. 

4. Modeling of Nozzle Flow Fields 

The analysis considered the effects of both changing flow area and heat 
exchanged. The effects of wall friction assumed to be negligible. The quasi 
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one-dimensional assumption will be used. Consider the steady flow through 
the control volume as shown in figure (2). 




Figure (2) One-Dimensional Control Volume 



Hence the continuity equation gives [10]: 
dp + dA + dV =Q ______ 

p A V 

While, the conservation of momentum gives [11]: 



P A + 



P + 



dP 



dA - (P + dP) (A + dA) = p A V dV 



(12) 



The second term on the left hand side represents the force due to the 
pressure on the wall. By neglecting the second order terms gives [10]: 

dP + pVdV = 0 (13) 

Also, the equation of state gives [10]: 

...... (H) 



dP dp dT 

= — - + 

P p T 



And the conservation of energy gives [11]: 

dq = C P dT (15) 

From equations (11) and (13) yields: 



dP-pV 2 



'dp dA A 

— ^ + 



v 



= 0 



P A , 

From equations (14) and (15) yields: 



(16) 



JO 

No.5 JOURNAL 

OF COLLEGE OF 
EDUCATION 2011 

dp _dP dq 



p P C P T 



(17) 



Substitute equation (17) into (16) and rearrange, leads to: 

'dA dq' 

dP = pV 2 ^ — %^ (18) 

p 

This equation gives the variation of pressure across the control volume, it 
mean that: 

P 2 =P,+dP (19) 

Where, 1 :- refers to the inlet of the control volume 

2:- refers to the outlet of the control volume 
From equation (15) get: 

T 2 =T 1+ ^_ (20) 

C p 

Use simple equation of state for perfect gas to determination the density 
value, hence: 

p 2 =-^- (21) 

R T 2 

The velocity value then can be found from the continuity equation, so: 

o 

V 2 = -5- ------ (22) 

P 2 A 2 

Where, A 2 = A x + dA (23) 

m =p 1 V 1 A 1 (24) 



The value of heat added that calculated in equation (9) is in (Watt) so it must 
be divide on the mass flow rate ( m ) to satisfy the requirements of the nozzle 
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modeling which is in (J/Kg). The output power then can be calculated at the 
exit of the nozzle by: 

P out =0.5pAV 3 (25) 

Demonstrating the critical length of the nozzle need to achieve the 
conversion of full solar energy absorbed into the kinetic energy of flow at 
the throat of the nozzle, it means: 

I s A N =0.5pA t V t 3 (26) 

Then the critical length can be determined as: 

L = p A ' V ' (27) 

&t *I,(D i+ D e ) U7) 

For air at operation design, the critical length is calculated to be (8 m); 

therefore the design length (10 m) is acceptable. 

5. The Mathematical Procedure 

In order to calculate the flow properties inside the nozzle, an program of 
FORTRAN-90 is adopted and explained in flow chart as shown in figure (3). 
The nozzle is divided into 40 control volumes, the boundary conditions (dq) 
is provided in each control volume and the inlet properties are taken as 
shown in table (1). The nozzle is made from steel, the inlet diameter is (10 
m), the exit diameter is (1 m) and its length is (10 m), and locates in 
horizontal level (7 m) up of the ground. Equations (18 - 24) are used to 
calculate the properties in the exit of the control volume then take the next 
control volume and repeat this procedure until reaching to the end of nozzle. 
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Cal. qr - Eq. (1) & Cal. ql - Eq. (4) 



Calculate q - Eq. (9) 



(N-l))?hdq = q/( 



5 




Calculate dp Eq. (18) 



Calculate P2 Eq. (19) 
Calculate T2 Eq. (20) 
Calculate R2 Eq. (21) 
Calculate V2 Eq. (22) 
Calculate A2 Eq. (23) 



i 
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Table (1) Local Air Properties 



Velocity 


3 


c p 


1005 




m/s 




J/Kg. 




ec 




K 


Pressure 


100 




2x 10" 




KPa 




5 








Kg/m. 








sec 


Temperature 


35 


K 


0.0268 




°C 




W/m. 








K 


Density 


1.17 


Pr 


0.7 




Kg/ 








m 3 







6. Results and Discussion 

Analysis of the design conditions and the interference effects of the free 
stream velocity and the heat load have been done. The flow velocity and 
temperature are the dominated parameters for comparison at different 
conditions. For confidence of understanding, figure (4) shows the variation 
of flow velocity along the nozzle, while figure (5) shows the variation of 
flow temperature along the nozzle. 

From figure (4), it is observed that the velocity of the air increase by the 
decrement of the area. The raise in the value of velocity in the case of heat 
added is reached to (50 m/sec) and in the case of no heat it is about (48 
m/sec), while Williams [6] determined that it is about (52 m/sec) for case of 
vertical nozzle. This mean that at heat gained the power reached to (57 KW), 
while in the case of no heat it is (54 KW). 

From figure (5), it is observed that the temperature of the air increased 
rapidly in the case of heat added, while it may be less in that of no heat. The 
temperature difference in this study is about (0.9 °C), while Williams [6] 
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determined that it is about (1.3 °C) for case of vertical nozzle and solar 
absorber arranged in the base of the nozzle. 

60.00 1 




o.oo — I 1 1 1 1 1 1 1 1 1 1 

0.00 2.00 4.00 6.00 8.00 10.00 

Distance From Start (m) 



Figure (4) Variation of Flow Velocity along the Nozzle 

309.60 1 




307.60 1 1 1 1 1 1 1 1 1 1 1 
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Figure (5) Variation of Flow Temperature along the Nozzle 



No.5 JOURNAL 

OF COLLEGE OF 
EDUCATION 2011 

Figure (6) show another comparison, it is observed that the percentage error 
of flow temperature between the experimental work of reference [7] and this 
study is about (-2 %) when (x < 7 m) and about (+4 %) when (x > 7 m). 

40.0 1 




Distance from Start (m) 

Figure (6) Comparison of Flow Temperature with Ref. [7] 

Although that the heat gained is a small, but it can be increase by advantages 
of concentrating the solar radiation [2]. The concentrating of radiation can 
be rise the solar power to about (4500 W/m 2 ), that's mean that the exit 
velocity in the design operation conditions reaches to an increment of (20%) 
than that of no concentrating (450 W/m 2 ). 



7. Conclusions 

Generally, there are several results to be noticed from the whole research 
according to the fluid energy which can be listed below as: 

1- It is an open cycle. 

2- A prototype solar nozzle is described of 10 m height, 10 m base 
diameter and 1 m throat diameter. 
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3- The value of velocity in the case of heat added is reached to (50 m/sec), 
while in the case of no heat transferred it is about (48 m/sec). 

4- Maximum heat gained could give an output power reached to (57 KW). 

5- The temperature difference through the nozzle is (0.9 °C). 

Furthermore analysis can be made by using of solar absorber in several 
layers thick arranged around the nozzle. Also put an attention to make a 
closed cycle with relatively by-passing, due to the difficulty of heating the 
air and an economic procedure. 
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